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We reconsider the calculation of 0(Is^ CD /ml) nonperturbative corrections to 
B — > X s l + l~ decay. Our analysis confirms the results of Ali et al. for the dilepton 
invariant mass spectrum, which were in disagreement with an earlier publication, 
and for the lepton forward-backward asymmetry. We also give expressions for 
the 0(k 2 Q CD /ml) corrections to the left-right asymmetry. In addition we discuss 
the breakdown of the heavy quark expansion near the point of maximal dilepton 
invariant mass q 2 and consider a model independent approach to this region 
using heavy hadron chiral perturbation theory. The modes B — > Kl + l~ and 
B — > KttI + 1~, which determine the endpoint region of the inclusive decay, are 
analyzed within this framework. An interpolation is suggested between the region 
of moderately high q 2 , where the heavy quark expansion is still valid, and the 
vicinity of the endpoint described by chiral perturbation theory. We also comment 
on further nonperturbative effects in B — > X s l + l~ . 



1 Introduction 



The inclusive decay B — > X s l + l~ (I = e, //, r) has received considerable interest 
in the literature As a loop-induced flavour-changing neutral current 

(FCNC) process it provides a sensitive probe of flavour dynamics, the least tested 
sector of the Standard Model. The rare decay modes B — > X s l + l~ are well within 
reach of the next generation of precision B physics experiments and promise to 
yield much needed information complementary to that from other sources such as 
B — > X s 7, B — > X s vu, B — B mixing, CP violation or rare K decays. The interest 
in B — > X s l + l~ and other inclusive rare B decay processes is reinforced by the fact 
that their theoretical treatment is fairly well under control. Indeed, the rate for 
B — ► X s l + l~ is dominated, in the region of q 2 = (pi- +P1+) 2 away from resonance 
backgrounds, by perturbatively calculable contributions. These are known at 
next-to-leading order (NLO) |7|, |[. Such a calculation at the parton level is 
formally justified by the heavy quark expansion (HQE), in which the free b quark 
decay b — > sl + l~ emerges as the leading contribution to B — > X s l + l~ . This result 
receives power corrections of the form {h.QCD/'m>b) n , which can be systematically 
addressed within the HQE framework. The leading corrections arise at order 
n = 2. They have been first considered in H. A further computation of these 



effects in |14| did not confirm the results obtained in ||. In particular, in JT4 
the relative 0(Aq CD /ml) correction diverges at the high-g 2 endpoint, indicating 
a manifest breakdown of the heavy quark expansion, a feature that is absent in 
§■ 

The C(Ag C . D /m 2 ) effects are relevant both conceptually, for assessing the validity 
of the HQE, as well as for obtaining quantitative control over a class of theoretical 
uncertainties beyond perturbation theory. In view of this, the phenomenological 
interest of B — > X s l + l~ and the situation in the literature described above, a 
further independent analysis of the issue is certainly useful. The results of such 
an analysis of the 0{A 2 ^ CD /m 2 ) corrections in B — > X s l + l~ will be presented 
in this paper. We will furthermore discuss the breakdown of the HQE near the 
endpoint of the spectrum (at maximum q 2 ). The implications of this feature for 
a description of the high-g 2 region will be pointed out. A major part of this work 
will then be devoted to investigating the model independent constraints on the 
g 2 -spectrum that can be obtained using heavy hadron chiral perturbation theory 
(HHChPT). 

The paper is organized as follows. A brief description of the general framework 
and a collection of basic formulas is given in Section 2. In Section 3 we present 
our results for the 1/m 2 corrections to the dilepton invariant mass spectrum, 
the forward-backward (FB) asymmetry and the left-right (LR) asymmetry in 
B — > X s l + l~ . Section 4 contains a discussion of the breakdown of the HQE near 
the endpoint. In this section we also analyze the endpoint region of B — > X s l + l~ 
in terms of the exclusive modes B — ► Kl + l~ and B — ► Kirl + l~ , calculated within 
chiral perturbation theory. A few comments on further nonperturbative effects 
in B — > X s l + l~ are made in Section 5. We summarize our results in Section 6. 
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2 Framework and Basic Expressions 



The starting point for the analysis of B — > X s l + l is the effective Hamiltonian, 
in the Standard Model given by (neglecting the small contribution ~ V* s V u b) 



G F 



8 a ~ - a ~ - 

Y^Ci^Qi + —C 9 (fi)(sb) v _ A (U)v + —C 10 (sb) V - A (ll) A 

.1=1 



(1) 

The Hamiltonian is known at next-to-leading order |7|, A detailed review may 
be found in [f2(J , where the Wilson coefficients and the four-quark operators Qi 
are defined explicitly (the operators are typically of the form Qi ~ (s6)(cc), for i = 
1, . . . , 6, whereas Q 7 ~ em b sa^ v {l + 75)6^ and Q s ~ gm b sa^ v {l + 7 5 )A a &G^J. 
From ([!]) the following general expression can be derived for the differential decay 
rate 

dT(B -> X.M-) _ G 2 F m\ a 2 3 m b 

dxdyds 192vr 3 ' ts tbl 4n 2 4nm 2 M B { ' 

L% \(\C e 9 ff \ 2 + \C 10 \ 2 ) Wf + 4m 2 \C 7 \ 2 Wr + 4m b Re C 7 C e 9 ff * W&} 



+L% {2Re Cl ff *C 10 Wf + 4m b Re C 7 C* 10 Wtf} 



Here m b (Mb) is the 6-quark (B meson) mass. Cg is a (scheme invariant) effec- 
tive Wilson coefficient that includes, in addition to Cg from ([!]), the contributions 
from the b — > sl + l~ transition matrix elements of 4-quark operators Q\, . . ., Qq. 
Next 

L % = PxuPiv + P2tf>iu ~ g^vPi ■ P2 and = -ie^pfp^ (3) 

are the symmetric and antisymmetric leptonic tensors, respectively (pi (p 2 ) is the 
momentum of l~ (l + ) and e 0123 = +1). We also set s = q 2 /ml (q = p x + p 2 ), 
x = 2p ■ pi/m 2 and y — 2p ■ P2/171 2 , where p is the 6-quark momentum defined 
as p M = m b v M in terms of the 5-meson four- velocity t> M = p b /Mb- The hadronic 
tensors W^ u can be written as W- IU = 21m T^ v where 

TV = 1 Jd'xe-^(B\Tjr(x)f 9 (0)\B), (4) 
T 97 = * J d*xe-^(B\Tjr(x)tf(0)\By-£, (5) 



fx e-^(B\T j^(x)j^(0)\B)^ 

q 2 

Tr = i [d 4 xe-^(B\Tj^(x)jr(0)\B)^, (6) 
J q 

jS = SY(l - 75 )fc , #" = sa^(l + l5 )b . (7) 

Here the B meson state \B) is taken in conventional relativistic normalization 
(B\B) = 2EV (the explicit appearance of the factor 1/M# in (|2]) is due to this 
definition). 

Evaluating the hadronic tensors to leading order in the heavy quark expansion, 



2 



eq. (^) reproduces the well known quark-level results for the B 
distributions and asymmetries. For instance, defining 



X s l + l decay 



R{s) 



£t(b -» x s m- 

T(B -> X c ev) 



one obtains upon integrating over x and y 



R(s) 



a 

471-2 



Vt 



t.s 



V, 



cb 



f{z)K{z) 



(9) 



X 



2s) 



/W7|2 



+ |C; 



10 



+ 4 



2 

1 + - 

s 



|C 7 | 2 + 12C 7 ReC' 9 e// 



Here f(z) = 1 — 8z 2 + 8z 6 — z 8 — 24z 4 lnz is the phase space factor and 
the QCD correction factor (z = m c jm^) entering T(B —>■ X c ev)\ k(z) can be 
found in |20|. Note that for the dilepton invariant mass spectrum R(s) only the 
symmetric part in (0) (proportional to L s ) contributes. In (PI) we have neglected 



flU/ 

0(mf/ml) and 0(m 2 s /ml) terms, as we shall do throughout this paper, unless 
stated otherwise. The expressions given are therefore applicable to the cases 
I = e, fi. The extensions of (|9]) to the case mj ^ (relevant for I = r) and 
m s are given in ||, |l2j. Neglecting the strange quark mass is a very good 
approximation except near the q 2 endpoint. This region, however, suffers from 
large nonperturbative corrections and the entire partonic approach has to be 
reconsidered there (we will come back later to this point). 

A quantity closely related to R(s) is the left-right (LR) asymmetry, which 
measures the difference in the rates of producing left handed or right handed 
leptons in B — * X s l + l~ decay. As discussed in [21, 
be directly extracted from (01). Defining 



22], the LR asymmetry can 



R L ' R {s) = R{s) 



°9 "~ * 2 



one has 
A L r{s) 



R L (s) 



R R 



a 

An 2 



ts 



V 



cb 



8) 

f( Z )K(z) 



(10) 



:n) 



;i + 2s) (-2(7i Re C!/ f ) - 12C 7 C W 



Another interesting observable that can be studied in B — > X s l + l decays is the 
forward-backward (FB) lepton asymmetry M, which can be defined as 

d 2 T(B -> XJ+l- 



A FB (s) 



1 

T(B -»• X r ev) 



dcos9 



ds dcos8 



-sgn(cos( 



(12) 



where 9 is the angle between Z + and B momenta in the dilepton center-of-mass 
frame. As shown in |TJJ Afb{s) is identical to the energy asymmetry introduced 
in 



Till . The NLO perturbative result for Afb(s) is given by 

1-s) 2 



Afb(s) 



3a 2 V t , 2 



47T 2 



V 



cb 



f(z)n(z 



-Re{C* l0 



2C 7 + s C e 9 ff ] } 



(13) 
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Interestingly enough, both Alr(s) and Afb(s) are sensitive to the relative signs 
between CV, Cg and Ci . These asymmetries therefore offer useful additional 
information on the underlying short distance physics. 



G(A 2 



A 



QCD 



mi 



Power Corrections to R, Alr and 



FB 



The hadronic tensors W^ v in (0) can be systematically expanded in inverse powers 
of the heavy quark mass using the operator product expansion (HQE) approach 
supplemented by heavy quark effective theory (HQET). The general procedure 
is described in great detail in [^J for the case of B — > X Ufi lv decay. The first 
corrections to the parton result (0(1)) appear at C(Ag CD /m^). To this order we 
obtain the following expressions for the hadronic tensors (after contracting with 

) 



ATrrribMB 



4vrM, 



B 



J dxdy L%Wr 



J dxdy L%Wg 



2ml 



[l-sY(l + 2s) 



3A 2 



2m\ 



[1 - 15s 2 + 10s 3 ) 



1 + 



2ml 



(14) 



Here 



Ai 



3m b 
4vrM ; 



B 



J dxdy L%Wr 



(B\h(iD) 2 h\B) 



2M 



A, 



B 



2m\ 



(5 + 6s - 7s 2 



1 + 



XT 
2ml 



:i - s) 



3A 2 6 + 3s - 5s 3 
2ml s 



1 (B\hga ■ Gh\B) 
6 2M 



-Ml 



B 



(15) 



(16) 



(17) 



with h the 6-quark field in HQET. 

The results in (|l^)-(|T6"|) agree with [ TJ] but differ from the findings of || . The 
contribution involving Wif is the same that appears in the case of semileptonic 
B -> XJu decay. Integration of (|TJ) over s yields (1/2) [1 + (Ai - 9A 2 )/(2m 2 )], 
reproducing the well known correction factor derived in p4fl . Inserting ([14]) - 
(|T6D into (TJ) we obtain for the 1/m 2 corrections to R(s) in (Q) 



3A2 


' a 2 


V u 


2 1 r 


2m\ 


^47T 2 


v cb 





'l - 15s 2 + lfV 



IC , 9 // | 2 + |C'io| 2 ) 



(6 + 3s - 5s 3 ) 



4|C 7 | 2 



(5 + 6s-7s 2 )4C 7 ReC' 9 e// 



+ W) R{S) 



;is) 
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Here we have used the normalizing semileptonic rate including terms of order 
l/ml 



T(B -> X c ev) 



1927T 3 



V cb \ 2 f{ Z )K{z) 



Ai 3A 2 g(z) 



2m 2 2m 2 f(z) 



g (z) = 3 - 8z 2 + 24z 4 - 24z 6 + 5z 8 + 24z 4 In z 



(19) 



(20) 



that can be found for instance in [^3 ] . Note that the correction due to the kinetic 
energy of the 6-quark ~ Ai is given as a simple overall factor (1 + Ai/(2m 2 )) for 
both B — > X c ev and B —>■ X s l + l~ and therefore drops out in the ratio R(s). 
Since, in contrast to A2, the quantity Ai is not well known anyway, its absence in 
flTBp is a welcome feature. 

Given the results in (p!4|)-(p^) it is straightforward to write down the 1/m 2 cor ~ 
rection for A LR (s) in (|TTf) 



$i/ m 2 A LR (s) 



3A 2 / a 2 
2mf 



ts 



cb 



f(z)K( Z ) 



+ (5 + 6s - 7s 2 )4C 7 C 



10 



+ 



1 - 15s 2 + 10s 3 ) (-2Cio ReC* ff ) 



■A L r(s) 



(21) 



This correction has been discussed previously in pl| , however based on the in- 
correct results of 0. 

To calculate the FB asymmetry it is necessary to contract Wg I/ and Wgj with 
the asymmetric component of the leptonic tensor. The relevant terms, expanded 
up to 0(A 2 QCD /ml), are given by 



2nm h 



— — - I dxdy sgn(w — x) L^.Wq 1 ' 
i b M B J y 6 vy 1 v 9 



s(l-s) 2 + -^s(3 + 2s + 3s 2 
6m^ 



A s 



2nM. 
leading to 

$l/m 2 b A FB (s) 



1 f 

— J dxdy sgn(y - x) L^ U W< 



2m 2 h 



s(9 + 14s - 15s 2 ) , (22) 



97 



Ai 



;i-s) 2 + -^(3 + 2s + 3s 2 ) 
6m^ 



A, 



2m 2 h 



(7+ 10s -9s 2 



(23) 



3Ao / a 2 



2ml \4tt 2 



V, 



V, 



cb 



+ (7+10s-9s 2 )2C 7 



+ 



f(z)K( Z ) 

g(z 



RelC 



10 



A F b(s) + 



s(9 + 14s - 15s 2 )C 9 
4Ai s 



eff 



3ml (1 — s) 2 



A FB (s) . (24) 



Also in this case our finding is in agreement with ||14|| . 
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4 The High-g 2 Region 



4.1 Generalities 

The size of the 0(A 2 Q CD /ml) corrections in flUf), ( f2lf ) and fl24]) is quite moderate 



at the level of several percent, for values of s below about 0.6. On the other hand, 
when s approaches the endpoint (s = 1), the corrections for R, A LR and A FB tend 
towards a nonzero value, while the leading term of these quantities vanishes as 
(1 — s) 2 . The relative correction thus diverges in the limit s — > 1 and the rate R(s) 
becomes even negative for s close enough to the endpoint. Obviously, the HQE 
breaks down in the endpoint region. From the expressions given above one may 
recognize that in the case of R(s) and Am this behaviour is exclusively related to 
the A2-term, whereas in the case of Apb also the kinetic energy correction ~ Ai 
is not well behaved in the limit s — > 1. We remark that these features are not 
shared by the result given in |J and have been first observed by the authors of 



14 



In order to account for nonperturbative effects that elude the HQE approach, 
fl4}| supplement the partonic calculation with a Fermi-motion model to predict 
the q 2 spectrum and the shape of the FB asymmetry in the entire physical region 
including the endpoint. Although such an approach could be useful in principle, 
in particular when employed in conjunction with experimental data (used e.g. to 
fit model parameters), we will not perform such an analysis here. Instead, we 
would like to discuss to what extent model independent predictions can be made 
for the B — > X s l + l~ spectrum. We will thereby focus our attention on the high- 
q 2 region. For this purpose we shall first discuss the nature of the breakdown of 
HQE in slightly more detail. We follow here the general discussion presented in 

mm- 

The central element in the operator product expansion of the tensors T^ v 
in (f|) is the s-quark propagator. This propagator emerges in the evaluation 
of the time-ordered products (Q) and determines essential features of the l/m& 
expansion. Denoting k = m^v — q, = — igA^ and neglecting the s-quark 
mass, the s-quark propagator in a gluon background field may be written as 

Ss{k) = k 2 + 2ik- D- ip + ie ' (25) 

Up to terms of order Aqcd = A, k is the momentum of the final state hadronic 
system. In the usual case, that is away from singular kinematical points, one has 
k ~ 777.6, k 2 ~ m 2 and consequently the hierarchy k 2 ~ m 2 . ^> k ■ D ~ nib A ^> 
p p ~ A 2 . Therefore one can expand S s (k) = fc/k 2 + 0{A/m^) and the usual 
HQE is valid. 

A different situation arises in the endpoint region of the lepton energy spectrum 
in B —>■ X CtU lv and for the photon energy spectrum in B — > X s ^(. Here one still 
has k ~ nib in terms of components, however the kinematics is now such that 
k 2 ~ m-feA. For the quantities in the denominator of Q25| ) this implies k 2 ~ m^A w 
k ■ D ~ rribA ^>JD p ~ A 2 . An expansion in A /nib is still possible, but k ■ D/k 2 
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is now of 0(1) and the corresponding effects have to be resummed to all orders. 



This is the case discussed in detail in p6| , [27j] for B — > X C)U lu and in |2§] for 
B -> X s7 (see also ||). 

We would like to stress that the situation encountered in the endpoint region of 
the q 2 spectrum in B — > X s l + l~ is substantially different from the two cases just 
described. For the kinematics that is relevant here, q 2 ~ m 2 ~ M|, it follows 
that k ~ A and k 2 ~ A 2 . Then all three terms in the denominator of (|25|) 



are of the same order of magnitude ~ A 2 . The heavy quark expansion breaks 
down completely and not even an all-orders resummation, of the type useful for 
B — > X C)U lu and B — > X s 7, can be performed. This conclusion is clear on physical 
grounds, since at q 2 ~ M\ the two leptons are emerging back-to-back, carrying 
almost all the energy released in the decay of the B meson. The final state 
hadronic system has very low momentum and we are in a regime of manifestly 
nonperturbative QCD. 

At this point we would like to emphasize a conceptual consequence of this 
discussion for the treatment of the q 2 spectrum in B — > X s l + l~ within a Fermi- 



motion model, as employed in In the case of the photon energy spectrum 
in B — > Xg'-f or the lepton energy spectrum in B — ► X c ^ u lv the resummation of 
leading singular contributions in the HQE leads to a description of the endpoint 
region in terms of a shape function [^7], [28[ [25]. The shape function depends 
on nonperturbative physics that can qualitatively, at least to some extent, be 
modeled by a Gaussian description of the 6-quark momentum distribution inside 
the B meson (Fermi-motion). As explained above, a similar interpretation does 
not exist in the case of B — > X s l + l~ . A Fermi-motion description of nonpertur- 
bative effects, particularly for high q 2 , appears therefore certainly less justified 
than in the usual applications to B — > X^ and B — > Xlv. In fact, as we have 
seen above, the divergence of the 1/m 2 corrections to the q 2 spectrum near the 
endpoint arises from the chromomagnetic interaction term ~ A2 that does not 
have an obvious interpretation in terms of a Fermi-motion ansatz. 

On the other hand, the kinematical situation in B — » X s l + l~ near the q 2 
endpoint, with few, low-momentum hadrons in the final state, lends itself to a 



treatment using heavy hadron chiral perturbation theory (HHChPT) 30 ] . 
Combining this description at very high q 2 with the standard HQE result at 
somewhat lower q 2 , where the latter is still valid, a model independent analysis 
of the entire high-g 2 region (above the \1/ and resonances) could be conceived. 
In the following we shall examine such a possibility. 

First, one may write down an effective Hamiltonian, suitable for the endpoint 
region (q 2 M^) in B — > X s l + l~ . This Hamiltonian differs from the standard 
Hamiltonian ([[]). 'Light' quark (u, d, s, c) loops may be integrated out explicitly 
since they involve the hard external scale q 2 ~ m 2 ^> lGeV. This endpoint 
effective Hamiltonian then takes the form, valid at NLO in QCD 

n eff , EP = -^V*V tb - x (26) 
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C 9 M^)v-a(11)v + C 10 (sb) v „ A (U) A + 2m b C 7 sa^(l + lh )b % M lv l 

The Wilson coefficient Cg^p has the structure 

C 9 EP = C$ DR + h(z, s)(3C[ 0) + C 2 (0) ) + (penguin contributions) , (27) 
with C[°\ C%\ Cg DR from (|l|). These quantities, the function h(z,s) and the 



remaining terms can be found in pPfl . Cg ^p is identical to Cg in (|^) (see also 



201), except that it does not include the QCD correction r](s) to the matrix 
element of the current (sb) v . A , which multiplies C 9 NDR in Cg". 

The Hamiltonian (^) is still normalized at a scale [i = 0(rrib). A further evo- 
lution down to hadronic scales ~ 1 GeV is calculable perturbatively using HQET 
('hybrid renormalization'). However the HQET logarithms will be automatically 
contained in the matrix elements of the {sTb) operators if they are taken in full 
QCD and appropriately normalized at fi = 0(rrib). It is therefore not necessary 



to make these effects explicit in fl26|). 



4.2 B Kl+l- 

At the very high end of the spectrum, between the Ktt threshold and the physical 
endpoint, the inclusive decay B — > X s l + l~ degenerates into the exclusive mode 
B — > Kl + l~ . Introducing the variable s m = q 2 /M B , this region corresponds to 
< s m < s* where S* = s m , max = (1 -M K /M B f = 0.821 and s** = 0.774. 
The matrix elements needed for B — > Kl + l~ can be written as 

(K(p K )\s^(l- l5 )b\B(p)) =f4 q i)( p + PK ) fl + f_tf)(p-p K )^ , (28) 

(K(p K )\sa^b\B(p)) = -ia T {q 2 W K f - f K v») , (29) 

in terms of the form factors f±{q 2 ) and ar{q 2 )- The decay rate (normalized to 
the semileptonic width as in @) is then given by 

-±-B(B ^> Kl+l-) T (B d ) GlMl a 2 

* M * mJ " - X c ez/) " B(B-X c ei/) 192*3 |V ^ s| 4^ MSmj X 

x {4 (I^bpI 2 + IC'iol 2 ) + ^m 2 |C 7 | 2 - / + a T m b Re C 7 C* P ) , (30) 



M 2 

/i(s m ) = ((1 - Q+ s m ) 2 - As m f 2 , Q = t§ ■ ( 31 ) 



v 2 V~^' ■ «-^« y ■ 2 
where the phase space function f\ reads 



In general, the form factors are very difficult to calculate. However, as long as 
we are interested in high q 2 , where the kaon momentum is small, HHChPT may 
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m c 




\Vts\ 


\v* 1 


4.8 GeV 


1.4 GeV 


0.04 


0.04 


1 


rh t {m t ) 


M w 


sin Bvk 


c^ 1 




170 GeV 


80.2 GeV 


0.23 


129 


0.225 GeV 


M B 


M K 


M n 


r(B d ) 


B(B -»• V c ez/) 


5.28 GeV 


0.496 GeV 


0.140 GeV 


1.6 ps 


0.104 


9 


f B 


/* 


A = M B * - M B 


//, = M Bs - M B 


0.5 


0.180 GeV 


0.132 GeV 


0.046 GeV 


0.090 GeV 



Table 1: Compilation of input parameters (central values). 

be used to estimate these nonperturbative quantities. In this approach, to the 
lowest order, one finds p9l BUI |32l |33|1 



/ fr(l±, Mb Z^ : ), v. tK = ^\f-^ (32) 

2/^ V w-^+A + ^y 2Mb 

U v-p K + A + n s 

Here /g and / w are the i? meson and the pion decay constants in the normalization 
where f n = 132 MeV. A = M B * — Mb = 46 MeV, /x, = M Bs — Mb = 90 MeV 
and g is the HHChPT parameter that determines B*Bn and D*Dir couplings at 
low energy. The value of g could in principle be inferred from a measurement 
of T(D* — > Dtt), but present data only allow to set an upper limit g exp < 0.7. 
According to the theoretical estimates of [[33] in the following we will assume 
0.4 < g < 0.6. 

4.3 B -► Krrl+l- 

Between Kn and Ktitt thresholds, i.e. for s^ 7 ™" = 0.728 < s m < s^ n , also the 
B — > K-kI + 1~ decay is kinematically allowed. No other modes are permitted and 
the hadronic invariant mass is still small enough to justify the use of HHChPT. 

The matrix element of the left-handed current relevant to B — > Kirl + l~ can 
be generally decomposed in terms of four independent form factors. Defining 

(K l (p K )vr J (p 7r )|s7 M (l-7 5 )6| J B(p)) = icy [ap^ + bp K ^ + cp^ - 2ih,£ m ^p 0l pP K pl K , 

(34) 

the lowest order HHChPT results are given by |34 

9f B M B 



fl v-p„ + A ' 

fB 



, (35) 



2g _ v-p n v-(p K -p n ) 



v ■ p n + A v ■ (p K + p n ) + n s 
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P K ■ P* - v ■ Pk v ■ 
[v -p n + A] [v ■ (p K + p n ) + ji s 



h 



9 2 h 



2/2 [ v . p „ + A)[v ■ (p K + p„) + A + fi a 



(36) 
(37) 



where |c_ 



c - 



2 1 coo I 



2\Cr 



1. 



We have checked the results (p5|)— (p7[) , first obtained by the authors of [34 



and agree with their findings. In addition we need the corresponding matrix 
element of the magnetic penguin operator. We obtain, again to leading order in 
HHChPT 



(K*(p K )^{p«)\sa^{l+ lh )b^\B(p)) 



(3f 



ic 



i'Pk,h + b'p^ + c'p^ - 2ih'e^ a ^p a p f3 K pl 



b' 



h! 



B 



M B - v ■ p K - v ■ p w 



f*q*(vp„ + A) 

V-p K V- (p K +p n ) ~P K -Pn 



+9- 



v ■ (p K + Ptt) + A + /i s 



g 2 f B M B p K ■ p n + Ml - V ■ Pit V ■ (p + Pn) 



ftf(vp n + A) 
gfB 

f%q*(v Pn + A) 



v ■ (p K + p n ) + A + fi s 
M B v ■ p-n - Ml - p K ■ p n 



(39) 
(40) 



+9 



Pk'PttV- (p K - Pit) - M 2 K v ■ p n + M 2 v ■ p h 



B 



V ■ (P K + Ptt) + A + fi„ 

M B -v -p K -v -pit 



2flq 2 (v- P7T + A) 



1 + 9 



v ■ (p K + p n ) + A + fi B 



(41) 
(42) 



We proceed to compute the decay rate. The necessary four-body phase space 
integrations can be performed using the general methods reviewed in P3 |. The 
leading behaviour of the differential B — > K-kI + 1~ decay rate as a function of 
(s^ — s m ) close to the Ktt threshold may be written down analytically. It gives 
the correct asymptotic behaviour at threshold and can be used as an approxima- 
tion to the full result for values of s m not too far from this point. We find 



Rk-tt+{St 



-^—B(B -> K-n+l+t 

ds m \ 

B(B -> X c ev) 



t(B a 



GlMl 



B(B -> Xreu) 192tt 3 



W a V u \ 



, a 
4^2 



■x 



X 



32tt 



l — {F 9 (s m ) (\C 9 , EP \ 2 + \C 10 \ 2 )+AF 7 (s m )\C 7 \ 2 + AF 97 (s m )Re C 7 C* BP } 



F 9 (s m ) 



,1/2 



IT (tlXxX2 
^(l-v^) 3 / 2 



< + 



1 " V^)^2 



(s Klt -s ) 3 



(43) 
(44) 
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with xi = x 2 = M K /M B , t x = (xi + x 2 ) 2 and s%* = (1 -Xi-x 2 ) 

functions F7 and F$ 7 are obtained from Fg by replacing wf 
respectively, where 



The 



wf and wf 



Wi 



f B M 1 



B 



J K 



M 



B 



M w + A \M K + M n 



1 + 



My + n s /2 
M K + M^ + fi t 



U ! 2 



__ f B M B gM n 



f B M B gM B 
' 2/2 Mw + A 

f B M B gM B 



(45) 



(46) 



w' 2 



fl M n + A M K + ' 

t 

Adding the two isospin channels, the total result for the nonresonant Bd(B~) 
KtcI + 1~ rate becomes 



2/2 M, • AM IS M K M., ■ 

(47) 



3 

Rk-k = Rk~tt+ + Rk°tt° — RrO-k- + RK-n° = tR-K~ 



(4? 



With the explicit expressions for Rk-k at hand, we are in a position to estimate the 
relative importance of the nonresonant Kit mode relative to the single K channel 
in the endpoint region. It is clear that the four-body process B — > Kirl + l~ is 
phase space suppressed. This is obvious from fl4"5D , which exhibits the typical 
factor of ~ 1 / (167T 2 ). More quantitatively we find that Rk-k amounts to less than 
2% of Rk at s m = 0.7 and is still less important for larger s m . Rk-k is therefore 
negligible in the entire endpoint region, which is completely dominated by Rk- 
The asymptotic formula (|43| ) describes the behaviour of Rk-k close to threshold 
(s m = 0.774). For s m = 0.7 (|43|) overestimates the full result by about 50%. This 
is still useful for an order of magnitude estimate. 

For s m below 0.7 a substantial enhancement of the Kir mode is expected due to 
the contribution of the K* resonance. However, for s m > 0.73 we are still far 
enough from the K* threshold to safely neglect the Kit mode with respect to the 
single kaon channel. 



4.4 Discussion 

In Fig. I] we compare the HQE result for R(s) with the HHChPT picture close 
to the endpoint. For this purpose we rescale the ratio R(s) (|S|) from quark to 
physical (hadron) kinematics, replacing 



R(s) - R(s m ) = -^R -^s m = . (49) 



d -B{B -> X s l + l 
m b V m b m J B(-B — > X c eu) 

This representation of the quark level result is furthermore useful since it makes 
the dependence of the prediction on the 6-quark mass explicit. The corresponding 
uncertainty, which will unavoidably exist in comparing theory with experiment, 
is illustrated in Fig. [I] for the representative range m b = (4.8 ± 0.1)GeV. In 
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0.5 0.6 0.7 0.8 

Figure 1: The dilepton invariant mass spectrum (dB(B — > X s l + l~) / ds m ) / 
B(B — > X c ev) = R(s m ) as a function of s m = q 2 jM\. For s m < 0.65 the 
NLO partonic calculation, including 1/m 2 effects, is used. There the lower, mid- 
dle and upper curves correspond to m&/GeV = 4.7, 4.8 and 4.9, respectively. 
For s m > 0.73 we show the HHChPT prediction for R(s m ), which is dominated 
by B — > Kl + l~ . Lower, middle and upper curve are obtained for g = 0.4, 0.5 
and 0.6. Linear interpolations between the two regions (0.65 < s m < 0.73) are 
indicated by dotted lines to guide the eye. The dashed curve illustrates a smooth 
interpolation using central parameter values. The thresholds for the various ex- 
clusive modes occur at s m = 0.821 (K), 0.774 (Kir), 0.728 (Kim), 0.691 ± 0.008 
(K*, ± half width). 

this context we recall that m;, here refers to the pole quark mass. In fact, since 
the NLO QCD calculation for b — > sl + l~ is available, the distinction of the pole 
mass from other mass definitions is already meaningful at first nontrivial (i.e. 
one-loop) order. The value of m& is to be determined from some other observable 
and can then be used as input for B — > X s l + l~ . In principle the error on 
can be further reduced in the future. We remark that the dependence of R(s m ) 
on the renormalization scale fi [m^jl < \l < nib) is less than ±5% in the region 
0.5 < s m < 0.7. 

The 1/ml corrections to B — > X s l + l~, which are included in Fig. J], are 
negative for s m > 0.5, increase with s m and reach about —20% of the leading 
result for s m = 0.65. As discussed above, nonperturbative effects that are beyond 
the control of the HQE become important for still larger values of s m . 

Very close to the endpoint at s m = 0.821 HHChPT offers a complementary 
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approach that may be used to constrain the behaviour of the spectrum from the 
region of large s m . An interpolation suggests itself between the regime s m < 0.65, 
where the HQE is valid, and s m > 0.73, where HHChPT may be used. In this 
way an essentially model independent description of the entire high-g 2 region 
s m > 0.5 (above the ^' resonance) could be obtained, at least in principle. In 
practice there are however several sizable uncertainties associated in particular 
with the HHChPT treatment. The B*Bit coupling g is still poorly known. Other 
uncertainties are related to the values of \V ts \ and r(B d ) entering (|3"0|), but these 
are less important than the one from g. Also the B meson decay constant fs 
introduces some uncertainty. 

A further issue is the reliability of chiral perturbation theory in the present case. 
The kaon mass is not very small with respect to the chiral symmetry breaking 
scale A x ~ 1.2 GeV. Thus, even in the vicinity of the endpoint, corrections of order 
30%-40% can be expected. In the Kti channel the situation could be even worse, 
given the presence of the nearby K* resonance. However, for a given value of s m 
the hadronic invariant mass ranges from to M™d = M#(l — s^ 2 ) ~ 770 MeV 
and only near the upper figure the effect of the resonance should be important. 
Given the above remarks, the result for R(s m ) we have presented should still 
provide a reasonable estimate. In addition, in view of the kinematical suppression 
of the Kn channel, the fact that the region above s m = 0.73 is entirely determined 
by B — > Kl + l~ can be expected to be valid beyond the limitations of chiral 
perturbation theory, which is useful for further studies. 

Systematic improvements are possible by going beyond the lowest order in 
HHChPT. In [32] chiral logarithmic corrections to the leading result have been 
investigated within HHChPT for the exclusive mode B — > Kl + l~. The correc- 
tions were found to be about 40%, which is sizable but still moderate enough for 
the approach to make sense. A related issue is the question of whether to use fx 
instead of which also goes beyond the leading order of chiral perturbation the- 
ory. The calculation of |I2| can be considered as a naive estimate of the expected 
size of the higher-order corrections, but lack of knowledge of the corresponding 
counterterms makes any precise statement about their exact value difficult. For 
this reason we have not explicitly included the chiral logarithms in our estimates. 
The related uncertainty is at least partly included in our variation of the coupling 
9- 

Apart from the differential branching fraction also the forward-backward a- 
symmetry can be studied in HHChPT at large q 2 . In this context we note that 
Aps vanishes identically for the single kaon mode B — > Kl + l~ . The endpoint of 
Afb is therefore determined by B — > Kiil + l~ and occurs at s m — (1 — (Mk + 
M n )/M B ) 2 = 0.774. 

We finally remark that the entire high-g 2 region (defined by 0.5 < s m < 0.821) 
corresponds to an integrated branching ratio for B — > X s l + l~ of about 0.5 ■ 10~ 6 
in the Standard Model. Thus, although the dilepton mass spectrum is dropping 
to zero towards the endpoint, a sizable branching fraction for B — > X s l + l~ exists 
in the region that is characterized by the transition from quark level dynamics to 
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HHChPT. The high-g 2 regime constitutes one of the interesting regions to search 
for B — > X s l + l~ in experiment |35], [36], [37]. Attempts to describe this part of 



the spectrum in a model independent way along the lines proposed in this paper 
should therefore be useful for the study of rare B decays at future B physics 
facilities. 



5 Other nonperturbative corrections 

In addition to higher order terms in the l/m& expansion, B —>■ X s l + l~ decays 
are affected by long-distance corrections related to cc intermediate states. These 
originate from the nonperturbative interactions of the cc pair in the process B — > 
X s cc — > X s l + l~ . If the dilepton invariant mass is close to one of the two narrow 
J PC = 1 cc-resonances (^(3097) and ^'(3686)) this effect is very large and 
"obscures" the short-distance FCNC process. However, this background can be 
eliminated by suitable cuts on the dilepton invariant mass. Given the vicinity 
of the two narrow resonances, two g 2 -regions naturally emerge as appropriate for 
the study of short-distance dynamics: the region below the ^ and the one above 
the \J r '. In the first case it is still necessary to deal with the cc rescattering below 
threshold, whereas in the second case the effect of broader resonances and open 
charm has to be evaluated. 

Nonperturbative contributions generated by cc intermediate states have been 
widely discussed in the literature by means of phenomenological resonance-ex- 
change models j|, [12|, [13], [TJJ. These approaches are useful near the main reso- 
nance peaks, but their validity outside this region is certainly less reliable. Indeed, 
the shape of the resonance tails far from the peaks is not under control. More- 
over, a double-counting problem is usually posed by the simultaneous use of quark 
and hadronic degrees of freedom. Within this framework, the only way to avoid 
double counting is represented by the approach of |12|] (KS). Here, in order to 
take into account charm rescattering, the correction to C$ induced by b — > ccs 
operators is estimated by means of experimental data on a(e + e~ — > cc-hadrons) 
using a dispersion relation. To be more specific, the function h(z, s) appearing in 
( p?D is replaced by 

h(z,s)^h(z,0) + -P ds'^f^ + i^R^s), (50) 
3 Js c s'(s' — s) 3 

where R^di 3 ) = o r (e + e~ — > cc)/a(e + e~ — > /i + /i~) and s c is the cc threshold. This 
method has also the advantage of including open charm contributions. However, 
it is exact only in the limit where the B — * X s cc transition can be factorized into 
the product of sh and cc color-singlet currents (i.e. non-factorizable effects are 
not included). Using this method we have estimated the long-distance corrections 
to the plot in Figure 1. The effect is quite small, at the level of several percent, 
essentially negligible for s m <; 0.53. Below this value the correction exceeds 10% 
because of the vicinity of the peak. 
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Larger effects from the higher cc resonances (\I/(3770), \I/(4040), \I/(4160), 
\I/(4415)) are obtained when a phenomenological factor k, pa 2.3 is introduced to 



enhance resonance production with respect to the factorization result [p8 |. This 



is motivated by the fact that the factorization assumption yields too small values 
for the B — > J/^>X S branching fraction. The validity of such a procedure for 
estimating the impact of higher resonances in B — > X s l + l~ is not entirely clear. 
Further work on this issue is necessary. In any case the deviations from quark- 
hadron duality due to resonances are reduced when the B — > X s l + l~ spectrum is 
integrated over a large enough range of q 2 . 

A more systematic and model-independent way to estimate cc long-distance 
effects far from the resonance region, based on a heavy quark expansion in inverse 
powers of the charm-quark mass, has been recently presented in fl6| (see also [|T7[| ). 
This approach, originally proposed in to evaluate similar effects in B — > 
decays, has the advantage of dealing only with partonic degrees of freedom. In 
this framework the leading nonperturbative corrections to R(s) turn out to be 
0(Ag CD /m^). They originate from the effective s6-photon-gluon vertex (induced 
by charm loops), where the gluon is soft and couples to the light cloud surrounding 
the b quark inside the B meson. The corresponding matrix elements can be related 
to A2 and thus are known both in magnitude and in sign. This kind of corrections 
is complementary to those computed in the KS approach, since they are generated 
by the charm rescattering in a color-octet state. Since the factorizable corrections 
vanish for s — > 0, as shown by ([50]), the 0(k 2 Q CD /m 2 ) effect is expected to be 
the dominant long-distance contribution for small values of the dilepton invariant 
mass. For s < 0.2 the relative magnitude of the correction is very small (at the one 
or two percent level). Higher-order terms become more and more important near 
the cc threshold, where the description in terms of partonic degrees of freedom is 
clearly inadequate. Using a simple order-of-magnitude estimate of higher-order 
terms, it has been shown that the leading corrections should provide a reasonable 
estimate of the effect up to s = ?>m 2 /ml ~ 0.26 (s m < 0.21 ) [[RJ. In this region 
the effect is below 4%. The 0(AQ CD /m 2 ) corrections are again very small above 
the \&' peak. 



6 Conclusions 

Within the framework of the heavy quark expansion we have computed the non- 
perturbative corrections of 0{K 2 ^ CD /m 2 ,) to the dilepton invariant mass spectrum 
and the lepton forward-backward asymmetry in B — > X s l + l~ decay. Our calcula- 
tions confirm the results of [Q for these quantities, which were at variance with 



earlier work ||. For completeness we have also written down the 0(AQ CD /mf) 
corrections for the lepton left-right asymmetry. 

In the main part of our paper we have then focussed on the region of high 
dilepton invariant mass q (with q 2 > Ml,). This is one of the relevant search 
regions in experiments looking for B — > X s l + l~ and corresponds to an integrated 
branching ratio of about 0.5 • 10~ 6 in the Standard Model. The HQE breaks 
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down for q 2 too close to its maximum value at the endpoint of the dilepton mass 
spectrum. This is signalled by a manifest divergence of the relative 0(h 2 Q CD /m 2 ) 
corrections in the limit q 2 — * m 2 , as already observed in [f[4 |. We have discussed 
conceptual aspects of this breakdown of the HQE for B — > X s l + l~ and emphasized 
that it is impossible to remedy the failure of the usual l/m^ expansion at the 
endpoint by an all-orders resummation, in contrast to the case of e.g. the photon 
energy spectrum in B —>■ X^. We were therefore led to consider an alternative, 
model independent approach to the endpoint region using HHChPT, which is 
in principle well suited in this kinematical regime. For this purpose we have 
formulated, at NLO in QCD, an effective Hamiltonian adapted to the endpoint 
region. This Hamiltonian is a variant of the standard Hamiltonian for b — > 
sl + l~ transitions and serves as the basis for calculating the relevant exclusive 
channels in the vicinity of q 2 = (Mb — M^) 2 within HHChPT. We explicitly 
considered the modes B — > Kl + l~ and B — > Kirl + l~ and demonstrated that the 
latter is completely negligible in the kinematical region of interest. To obtain 
a complete description of the high-g 2 spectrum, an interpolation between the 
HHChPT regime and the region of validity of the heavy quark expansion has been 
suggested. At present there are still limitations in accuracy from uncertainties in 
the value of rrib and, particularly, in the poorly known HHChPT parameter g as 
well as due to neglected higher order terms in the chiral expansion. However, the 
approach is essentially model independent and systematic improvements can in 
principle be made. 
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